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Abstract 

One can define several properties of wave equations that correspond to the absence of tails in 
their solutions, the most common one by far being Huygens' principle. Not all of these defi- 
nitions are equivalent, although they are sometimes assumed to be. We analyse this issue in 
detail for linear scalar waves, establishing some relationships between the various properties. 
Huygens' principle is almost always equivalent to the characteristic propagation property, and 
in two spacetime dimensions the latter is equivalent to the zeroth order progressing wave prop- 
agation property. Higher order progressing waves in general do have tails, and do not seem to 
admit a simple physical characterisation, but they are nevertheless useful because of their close 
association with exactly solvable two-dimensional equations. 
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1. Introduction 



The question whether a propagating wave leaves a tail behind it is one that has interested 
many physicists and mathematicians for a long time, and has found applications from the 
first studies in light propagation [1] to the theory behind the proposed experiments to de- 
tect gravitational waves [2,3]. Intuitively, a wave tail can be described as follows. Suppose 
that, at some time t, an instantaneous pulse of a field <p is produced at a point labelled 
by spatial coordinates {x 1 }. This pulse originates a wave front propagating outward with 
some speed, that will be detected by an observer sitting at some point {x n } at a time 
t' > t. We say that the wave has a tail or wake if such an observer continues to detect a 
nonvanishing field even after the wave front has passed, i.e., at times greater than t' . In 
this case we speak of diffusive, as opposed to sharp, propagation. 

One can roughly identify three distinct possible reasons for the occurrence of tails: 
(i) Mass-like terms in the wave equation. For example, the Klein-Gordon equation 



in four-dimensional Minkowski spacetime has tails when \i ^ 0, but not when \i = 0. 
This feature can be seen as the wave-mechanical counterpart of the fact that massive 
particles move slower than the speed of light. 

(ii) Dimensionality of spacetime. For example, the massless Klein-Gordon equation in 
m-dimensional Minkowski spacetime has tails for odd values of m, but not when m is 
even, with the exception of m = 2 [4-7]. 

(iii) Backscattering off potentials and/or spacetime curvature [3,8,9]. This is clearly the 
most interesting tail production mechanism from a physical viewpoint. 

The study of the occurrence of tails, as well as of their implications, is certainly of 
great interest and importance, but relies on the possibility of using precise definitions. The 
heuristic characterisation given above is evidently too loose for this purpose, but it can 
be easily formalised to obtain a tail-free property. For linear wave equations, however, 
a no-tails condition is better expressed in terms of the Green function, and corresponds 
to what is usually called the Huygens principle. Furthermore, still other definitions (e.g., 
characteristic propagation property, progressing waves propagation property) can be found 
in the literature. These characterisations are not all equivalent, although they are some- 
times (explicitly or implicitly) assumed to be. It is the purpose of this paper to explore 
this issue in detail, by establishing the relationships between them and clarifying their 
meaning and domain of applicability. 

If the wave equation is inhomogeneous (which corresponds, physically, to the presence 
of sources), the formulation of no-tails properties is a delicate task, because of the need to 
capture the notion of sharp propagation when the field due to the sources is superposed 
to the one arising purely from the boundary conditions. The situation is even worse for 
nonlinear waves. For these reasons, we shall restrict ourselves here to considering an 
arbitrary linear, second order homogeneous hyperbolic partial differential equation, which 
can always be written in the form 



H(f)-fj, 2 (f) = 
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^dn<t> + h»dJ + Kct> = 0, 



(1.2) 
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where g^ u , and K are suitable real functions of m variables that we can think of 

as coordinates over an m-dimensional differentiable manifold Ai; without loss of generality, 
we can assume that = g UfJ- so that these functions can be thought of as the components 
of the inverse of a Lorentzian metric g ab on M. given, in the coordinates {x^}, by the 
inverse of the matrix g^; this metric, however, may not be the one normally used to 
determine intervals and causal relations (for a physical example in which it is not, see Ref 
10). The metric g ab defines then a Riemannian connection V that allows us to rewrite 
(1.2) equivalent ly as 

g ab V a V b cf> + H a V a <f> + K ( f> = 0, (1.3) 

where H a is uniquely determined by /i M and g^ v (specifically, = d v g^ v — In y/—g, 
where g is the determinant of g^ v )- 

Although physically one may wish to discuss the appearance of tails in waves emitted 
by very distant sources, the phenomenon is essentially a local one. Therefore, it will not 
represent a loss in generality to study wave propagation close to the hypersurfaces where 
data are given, and we can assume that M. is globally hyperbolic and normal. A more 
general manifold can always be covered by a collection of regions with these properties, 
and tails will develop iff they already do inside some of those smaller domains [11]. In 
such a manifold, a very powerful tool for studying general properties of linear differential 
equations and their solutions is that of Green functions [6,12]. A Green function G(x,x') 
for (1.3) is defined by 

(g ab V a V b - H a V a + K- V a H a ) G(x, x') = -6(x, x') , (1.4) 

and by appropriate boundary conditions; in (1.4), b~(x, x') is the delta function on A4, such 
that for each test function /, 

J d m x' x/^tM f(x')S(x,x') = f(x) . (1.5) 

It should be noted that the differential operator acting on the first variable in G is the 
adjoint of the one which acts on <p in (1.3), and differs from it when H a ^ [6]. 

If we choose to study the behaviour of the field at points to the future of a Cauchy 
hypersurface S on which data are assigned, it is convenient to work with the advanced 
Green function, which satisfies G(x, x') = when x' is to the past of x. (Of course, the 
whole discussion can be carried out as well for the "time-reversed" situation.) In general, 
since the wave equation gives rise to causal propagation ([13]; see also Ref 14, p 250), the 
advanced Green function can be written as 

G{x, x') = t(x, x') + A(x, x') , (1.6) 

where S is a distribution term with support on light-like separated pairs (x,x'), and A 
has support on timelike related ones, with x to the past of x' in both cases (in fact, this is 
a possible definition of causal propagation). Thus, A is of the form 

A(x, x') =: A(x, x') 9A-T(x, x')) , (1.7) 
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where T(x, x') is the square of the proper distance calculated along the unique geodesic 
connecting x and x', and 9 , , with a common abuse of notation, represents an "advanced 
step function," which is zero when x lies in the future of x' [6,12]. 

We will allow the Cauchy hypersurface S to be piecewise smooth, like, e.g., the future 
null cone of a point, which is not differentiable at the vertex. The field at each x G 
D + {S) = J + (S) can then be expressed as 

<j>{x) = - J dS a (x') g ab (x')G(x',x)V^(p(x')+H a (x')G(x',x)(j)(x') , (1.8) 

where dS a (x') is the oriented volume element on the hypersurface S at x'\ this expression 
is valid even if S is a (partially) null hypersurface, provided that dS a is appropriately 
defined (see section 5). 

We start in section 2 with a general discussion of the two-dimensional case. In section 
3 we give some definitions of properties of wave propagation related to the absence of 
tails, whose meanings and mutual relationships are investigated in sections 4-7. Section 
8 contains some concluding remarks, as well as speculations about possible lines of future 
research on the topic. As we have already been doing, we will use latin indices a, 6, 
... as abstract indices in spacetime [14], which just indicate the tensorial nature of an 
object without requiring a set of coordinates, while greek indices fx, z/, ... will be used 
for equations valid in some chart. The notations D ± (A), J ± (A), I ± (^l), where A is 
some subset of M, stand, respectively, for its future/past domain of dependence, causal 
future/past, and chronological future/past [14], all defined in terms of the causal relations 
induced by g ab . Minkowskian coordinates, in which the coefficients of the metric have (at 
least at a point) the form r]^ v = jf" = diag(— 1, 1, . . . , 1), will be denoted by (t, x l ), with 
i = 1 m — 1. 



2. The two-dimensional wave equation 

We begin with a general discussion of the wave equation (1.3) in a two-dimensional space- 
time. The motivation for devoting an entire section to such a specific subject is twofold. 
First, it will provide us with a concrete example with which to illustrate some general 
ideas. Second, and more important, much of the paper will turn out to deal exclusively 
with this case; it is thus convenient to establish a few preliminary points, to which we can 
refer later on. 

Any two-dimensional spacetime is conformally flat, and its metric can therefore be 
locally written as g ah = Q 2 r] ab , where Q is a nonvanishing function and r\ ah the Minkowski 
metric. In a general m-dimensional conformally flat spacetime we have 

v a x a = d a x a + d a in\n\ m x a = \n\- m d a (\n\ m x a ) , (2.1) 

for an arbitrary vector field X a . In the particular case m = 2, and for X a = g ab V h f = 
Q~ 2 rj ab d h f ', equation (2.1) allows us to write 

g ab V a VJ = O- V 6 W , (2-2) 
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for any function /. Therefore, in 1 + 1 dimensions, 

9 ab V a V b f + H a VJ + Kf = IT 2 (v ab d a d b f + H a dJ + Kf) , (2.3) 

where H a := Q 2 H a and K := Q 2 K. Since Q ^ everywhere on Ai, we have that <p 
satisfies the wave equation (1.3) in (A4,g) iff it satisfies 

T] ab d a d b (t) + H a dJ + K(j> = (2.4) 

in the flat spacetime (A4, rj). We can therefore restrict ourselves to studying (2.4) without 
loss of generality. 

It is convenient to introduce locally on M. the null coordinates 

u:=\{t-x), (2.5) 

v:=±(t + x), (2.6) 

in which the Minkowski metric and its inverse have only one independent nonvanishing 
component each, r\ uv = —2 and r\ uv = —1/2, respectively, so that (2.4) becomes 

dlj + Ud u <p + Vd v( j ) + W ( j ) = 0, (2.7) 

with U := —H u , V := —H v , and W := —K. In the rest of the paper, we shall often refer 
to this form of the two-dimensional wave equation. 

It is not difficult to check, using (2.1) again, that for an arbitrary function / in two 
dimensions, 

<7 a6 V a V 6 / - V a (H a f) + Kf = -IT 2 [d 2 u J - d u (Uf) - d v (Vf) + Wf] . (2.8) 

Since S(u, v; u', v') = O -2 S(u— u') S(v—v'), it follows that G(u, v; u', v') is a Green function 
for (1.3) iff it satisfies the equation 

d 2 uv G - d u (UG) - d v (VG) +WG = 5{u - u') 8{v - v') . (2.9) 

By direct substitution, we can verify that 

G(u, v; u', v') = A (it, v; it', v') 9{u' — u) 6{v' — v) , (2-10) 

where 9 is the step function, satisfies (2.9) provided that the following conditions are 
fulfilled: 

d 2 uv A - d u (UA) - d v (VA) + WA = 0; (2.11) 

A(u, v'; u', v') = exp ^- J du" V(u", v')j ; (2.12) 

A(u', v; u', v') = exp ^- dv" U{u', v")j . (2.13) 

Since (2.12) and (2.13) can be regarded as data in a characteristic initial value problem for 
(2.11), which has a unique solution, it follows that A (it, v; u', v') is completely determined 
by (2.11)-(2.13); therefore, (2.10) is the general form of the advanced Green function for 
the two-dimensional wave equation. In the particular case U = V = K = 0, we recover 
the well-known result A = 1 [5,7]. 



3. Definitions 



Various properties are used to characterise the absence of wave tails. First of all, we need 
to distinguish between the fact that a particular solution of the wave equation may not 
have tails, and a possible intrinsically non-diffusive nature of the equation itself. We say 
that a set of Cauchy data with compact support C produces no tails iff the field 4> obtained 
evolving these data vanishes in K + (C), the set of all points in the causal future of C which 
cannot be reached from C by a null geodesic [15]. For such data and points x G K + (C), 
(1.8) becomes 

<j>{x) = -j dS a (x') [g ab (x') A(x', x)^cj>{x') + H a (x') A(x', x) cj>{x')] 9{-Y{x' ', x)) , (3.1) 

since the E term and the derivative of the ^-function in G(x' ', x) do not contribute to the 
field at x (we retain the factor 9(— T) in order to cover also cases in which x G K + {C) is 
not in the causal future of some x' G C, as it may happen, for example, when C is not 
simply connected). Therefore, the data on C produce no tails iff the integral on the right 
hand side of (3.1) vanishes. 

A related notion, expressed directly in terms of the field rather than of the Cauchy 
data, generalises the concept of distortionless propagation, without necessarily referring 
to a differential equation [5]. We say that a wave propagates without distortion if it can 
be written as 4>{x) = R(u(x)), where the phase u is a function on M., and the wave 
form R is a function of one variable. Clearly, if R has compact support [u 1 , u 2 ] , the wave 
propagates permanently sandwiched between those two values of u. We have here no 
metric with which to verify whether this wave has tails in the sense defined above, but it 
is nevertheless obvious that "there are no tails with respect to u." A slightly more general 
situation is that of relatively undistorted or simple progressing waves, those of the form 
4>{x) = f(x) R(u(x)), where the amplitude f is a function on M. [5,6,16-19]. An example in 
Minkowski spacetime is the spherical wave 4>{t, x) = expi(fc|x| — u;£)/|x|, which progresses 
with speed ui/k in the radial direction, and is undistorted except for the decrease in the 
amplitude l/|x|.t A finite sum of simple progressing waves will also be called (relatively) 
undistorted, the simplest example being the general solution (f)(t, x) = R(t — x) + S(t+x) of 
the two-dimensional equation (2.7) with U = V = W = 0. A useful further generalisation 
is that of progressing waves of order N, those that can be written as 

N 

t{x) = Y J fi(x)R {i) «x))> (3-2) 

where the superscript (i) denotes the i th derivative. For these waves it is still true that a 
compact support for R implies a "sandwich propagation" for 0. 



t Progressing waves become interesting for suitable choices of the phase; the most common ones in 
Minkowski spacetime have u(t, x) = t — F(x), for some wave front F, function of the spatial coordinates 

[16]. 
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Turning now to intrinsic properties of the wave equation, corresponding to the absence 
of tails in the general solution, the most frequently used one by far is Huygens' principled 
HP: A wave equation is said to satisfy the Huygens principle (HP) iff the field at a 
point x depends only on the data at the intersection of the past light cone through 
x with the Cauchy hypersurface [4-6], in the sense that any two sets of data which 
coincide there must give the same field at x. 

It is obvious from (1.8) that HP is equivalent to the requirement that the advanced Green 
function G(x, x') have support only on pairs of points such that x lies on the past light cone 
of x' . The latter statement is often also referred to as "Huygens' principle". An obvious 
consequence of (2.10) is that HP is always violated in a two-dimensional spacetime. 

The concepts discussed at the beginning of this section give rise, however, to other 
definitions that can be found in the literature. Using the notion of data that produce no 
tails, we define: 

TF: A wave equation is said to be tail-free (TF) iff each set of Cauchy data with 
compact support produces no tails [15]. 

It is almost trivial to see that TF and HP are equivalent; nevertheless, for the sake of 
clarity, we shall present the explicit proof in the next section. 

It is sometimes interesting to consider data assigned on a characteristic (i.e., null [6]) 
hypersurface, and to give therefore a modified version of TF. This property was originally 
formulated in terms of solutions of a wave equation [20], but we find it more meaningful 
to state it as applied to the equation itself [15,21]: 

CPP: A wave equation is said to satisfy the characteristic propagation property (CPP) 
iff each set of null data of compact support, which vanish at any singular points of 
the support, produces no tails. 

Why restrict, in this definition, the admissible data to those that vanish at singular points? 
For a two-dimensional characteristic initial value problem, we can consider without loss 
of generality the initial data <p{u, 0) =: ip(u) and 0(0, v) =: ip(v) to be assigned on the 
union of the two "half-axes" u > and v > 0. Consider the simplest example of (2.7), 
with U = V = W = 0. Its general solution is of the form 0(it, v) = R(u) + S(v), which in 
terms of the data becomes 0(w, v) = <p(u) + ip(v) — 0(0, 0) . This means that, if 0(0, 0) ^ 0, 
no data with compact supports [1*1,1*2] and [i>i,i> 2 ] can lead to "sandwich propagation," 
i.e., one for which the support of the full solution is contained in the union of the strips 
u 1 < u < u 2 and v 1 < v < v 2 - There is, of course, nothing bad in this, since we know 
already that HP is violated in two dimensions. However, data <p(u) and ip(v) such that 
(p(0) = -0(0) = 0(0,0) = do lead to sandwich propagation. The point (0,0) here is an 
example of singular point of S, and the reason we excluded data which do not vanish at 



T The improper use of the term "principle" to denote what is actually only a property is commonly 
adopted for historical reasons, with the probably unique exception of Hadamard, who refers to it as 
"Huygens' minor premise" [4]. In reality, Huygens made use of TF below rather than HP, and merely 
in an approximate version, by postulating that almost all the waves emitted by a pointlike source are 
concentrated on the wavefront (see Ref 1, e.g., pp. 18 and 22). It is interesting to notice that this 
hypothesis does not correspond to assuming that tails are absent, but only that they are small! 
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such points in the definition of CPP is that this allows us to give a reasonable no-tails 
characterisation that may hold even in cases in which HP does not. Actually, we shall 
see in section 5 that CPP, although in principle weaker than HP and TF, in practice 
always implies them, except in two spacetime dimensions. We wish to stress that the 
restriction of data we are considering is actually a very small one, concerning a subset of 
C of measure zero which contains points that are already pathological. Furthermore, such 
a restriction corresponds precisely to what is commonly done when assigning data on the 
asymptotic past (probably the only characteristic initial value problem with good physical 
motivations), where one allows 4> 7^ at past null infinity X~ but requires 4> = at past 
timelike infinity %~ [14]. 

In addition, the notion of progressing waves motivates the following definition for 
two-dimensional equations: 

PWjy : A wave equation in two spacetime dimensions is said to satisfy the progressing 

wave propagation property of order N (PW^) iff its general solution can be written 

as a sum of two progressing waves, 



where the amplitudes f i and g i are fixed functions on M. depending on the wave 
equation (and at least one of f N and g N is not identically zero), while the wave forms 
R and S are arbitrary sufficiently differentiable functions of one variable (this forces 
the coordinates u and v to be null) [17,19,22]. 

In section 6 we shall show that, in two spacetime dimensions, CPP is equivalent to PW ; 
as far as the PW^, with A > 0, are concerned, their resemblance to a no-tails property is 
only apparent. Furthermore, wave equations in more than two dimensions may well have 
progressing wave solutions for appropriate choices of the wave front (see, e.g., Ref 16), 
but it seems excessive to ask that their general solution be expressible as a finite sum of 
progressing waves. These issues, which diminish the appeal of PW^ as regards the study 
of tails, will be discussed in the concluding section. 

Another property that can be satisfied by two-dimensional wave equations is their 
solvability by the method of Kundt and Newman [20]. To apply this method, we start by 
writing (2.7) in the two equivalent normal forms 



where <f) Q and ip are related to <f> by factor transformations <p Q = (pexpa and ip = </>expr, 
with 



N N 



<K*) = E fi( x ) S^(v(x)) , 



(3.3) 



i=0 i=0 



0«UA0o) -3i </>o = > 



(3.4) 
(3.5) 




(3.6) 




(3.7) 
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while the coefficients are related to those in (2.7) by 



Jo = 1 = ex P( r - a ) » ( 3 - 8 ) 
j 1 = (d v V + UV-W) Jo , (3.9) 

and 

l^VJJ + UV-W)^. (3.10) 
If we inductively define j k and 4> k by 

h+i/Jk = k/k-i ~ d L ln \k\ > ( 3 - n ) 
jk+Ak+i = 3 k d u<t>k > ( 3 - 12 ) 

assuming of course j fc 7^ for all k G Z, we obtain from (3.4) a countable set of equations 

0fc) - 3k+i <Pk=Qi keZ. (3.13) 
Similarly, if we inductively define l k and ip k by 

= h/h+i ~ d L ln \h\ > ( 3 - 14 ) 
h-\^k-\ = h d v^ki ( 3 - 15 ) 

assuming now that Z fc 7^ for all fc G Z, we obtain from (3.5) a second countable set 

du(!kdv1>k)-lk-i1> k = Q, ^z. (3.16) 

We shall refer to equations (3.13) and (3.16) as being in the fc th i>- and it-normal form, 
respectively. It is not hard to check that for all k € Z the /c th w-normal form equation, 
corresponding to the coefficients j k and j k+1 , and the /c th u-normal form equation, corre- 
sponding to the coefficients l k and l k _u are equivalent under the transformation 

0* (3-18) 

for fc e Z. It is also easy to see that the equations within the set (3.13) (respectively, (3.16)) 
are equivalent in the sense that a solution of any one of them generates a solution of every 
other one of them through (3.11)-(3.13), (3.17), and (3.18) (respectively, (3.14)-(3.18)), 
and we thus obtain two equivalence classes of v- and w-normal form equations labelled by 
the index k ranging over both negative and positive integers. 

Given a wave equation in the th w-normal form (3.4), we say that its substitution 
sequence {j k } is double terminating in N steps [20] when j kl+1 = and l k2 _ 1 = (but 
j k 7^ and l k2 7^ 0) for some k 1 > 0, k 2 < 0, and N = max {k x , — k 2 }. However, it 
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was shown in Ref 20, using (3.11)-(3.18), that in this case the general solution of (3.4) is 
00 = 4>a + h<t>Ri where 



^4<HW-'¥4*- sw )---)))< (3 - i9) 
*«-£<b<£<"-^<'**->)---)))- (3 - 20) 

with S(v) and R(u) arbitrary functions of one variable. It is obvious from (3.19) and (3.20) 
that such a is a progressing wave of order N. This relationship motivates the following 
definition: 

KNjy-: A wave equation in two spacetime dimensions is said to be solvable by the 
Kundt- Newman method in N steps (KN N ) iff its substitution sequence is double ter- 
minating in iV steps. 

As we have just seen, all KN N wave equations are PW^. We shall see in section 7 that the 
converse is also true, so that KN N and PW N are equivalent properties. 

4. Equivalence between HP and TF 

We now begin studying the relationships between the various properties of wave equations 
we listed in the previous section. The first three are related in a simple way in any number 
of dimensions; as we will now show, TF and HP are always trivially equivalent, which 
justifies the fact that they are often identified, and they are almost always equivalent to 
CPP (see next section). These results generalise to arbitrary wave equations of the type 
(1.3) those of Ref 15. 

In order to show that TF =>- HP, let us consider any point x G D + (S), and assign data 
with support C C I~(x) fl S. Since TF holds by hypothesis, and x G K + (C), it follows 
that 4>(x) = for each set of data prescribed on such a C. By causality, <j>(x) cannot be 
influenced by data given outside I~(x) fl S, hence it can depend only on their value on 
E~(x) fl 5, where E~(x) stands for the set of points in J~(x) that are null related to x. 
This is precisely the content of HP. 

Let us now prove that HP TF. First of all, let us notice that, by choosing data 
that vanish everywhere on S, we get 4>(x) = everywhere in Ai by (1.3). If we now assign 
data on S with compact support C, HP implies that the value of (j)(x), for any x G K + (C), 
must be independent of the data; in particular, it must have the same value than in the 
case of vanishing data, i.e., it must be equal to zero, by the remark above. This completes 
the proof that HP ^ TF. 

It is interesting to notice that the structure of this proof allows one to extend it to 
more general wave equations than (1.3). Actually, in the case of the latter, a simpler proof 
can be given that makes use of the equivalence between HP and the property A(x', x) = 0. 
In fact, it follows immediately from (3.1) that, if HP is satisfied, then <f>(x) = for all 
x G K + (C), i.e., the wave equation is TF. The converse is also true, as we can see by 
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choosing S to be spacelike at a point x, and data of support only at x, such that, denoting 
by n a the unit vector normal to 5, 

cf>\ s (x) = 0, (4.1) 
(n a V a cf>)\ s (x) = 8 s (x,x), (4.2) 

where 5 S is the (m — 1) -dimensional delta function on S. Then (3.1) reduces to 

(j)(x) = -A(x,x) , (4.3) 

and if 4>(x) = for all x G K + (C) = I + (x), we must have A(x, x) = for all such pairs of 
points, i.e., A(x', x) = 0. However, this proof relies heavily on the existence of the integral 
representation (1.8), and appears much more difficult to generalise to the case of more 
complicated wave equations than the "geometrical" one given above. 



5. Relationship between HP and CPP 

It is obvious from the definitions given in section 3 that TF, and thus HP, imply CPP. To 
check whether the converse holds we need to consider again the integral expression (3.1) 
for x G K + (C), where C is now part of a piecewise null hypersurface S locally defined 
by the condition w = const, for some w : A4 — > R with g ab V a wV b w = at all points 
where S is differentiable (this condition may fail only in a subset of S of measure zero). 
One usually writes dS a = dS n a , where dS and n a are often thought of as volume element 
and unit normal to the surface <S, respectively, which clearly does not make sense if S is 
null. However, the definitions of these objects can be extended to that case as well. Let 
us introduce local coordinates {£* : i = 1, . . . , m — 1} on S, and use w as a coordinate 
transverse to S, so that (w,£ l ) are coordinates on Ai adapted to S. Then we can always 
write dS^ = d m_1 £ y/— g S^ w , with the understanding that g must be calculated in these 
coordinates, which in the null case can be split into dS = d m-1 £ yf^g, and = 5^; the 
latter are the components of the form n a = V a w. With these conventions, we can write 
the Gauss theorem for an arbitrary vector field Y a as 

/ dVV a Y a = I dS a Y% (5.1) 

J AT JdM 

where Af is a region in AA, and dV = d m Xy/—g is the spacetime volume element, even 
when (part) of dftf is null. This follows from the fact that the left hand side of (5.1) can 
be broken into a sum of integrals of the form 

/ d m -^dwd^(y/^Y^) , (5.2) 
Jv 
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where V is a domain in M m one of whose "faces" (that we denote by J 7 ) is defined by 
w = 0; but (5.2) can be transformed, using the Gauss theorem in M m , into a sum of terms 
corresponding to the various faces of V, of which only 



d^V^S^Y" (5.3) 



contributes, in the end, to the total expression. Therefore, (5.1) (hence (1.8) and (3.1)) is 
completely justified, independently of the causal character of dM . 
We can rewrite (3.1) in the coordinates (w,C) as 

m = - [ dS'{ V>"A0) - <j> [(V/)A + n"(2 V^A - H^A)] } 9(-T) (5.4) 

(here and in the following integrals, all functions in the integrand depend on x', and all 
two-point functions on the ordered pair (x',x)). Since we have = g^ w ', which implies 
n w = g ww = g^ v d w d v w = 0, the first term in the right hand side of (5.4) contains only 
derivatives performed along 5, and can be transformed as 



fd5'V;(/A^) = -f d m - 1 ed , i (V^n i Acj ) ) = -<f da> a A, 

JC ' JT JdC 



(5.5) 



where T stands again for the region of M m defined by w = 0, da a is the oriented volume 
element on dC, and the Gauss theorem in M m_1 has been applied in the last step (pos- 
sible discontinuities of <fi do not prevent us from using the theorem, if they are correctly 
interpreted in the sense of distributions). We can therefore deduce that this term does 
not contribute to the expression, because <f>(x') = outside C and we can think of the 
integration as being performed over a region larger than C, on whose boundary <f>(x') = 0. 
Therefore, CPP is satisfied iff 



dS' <P [(VX) A + (2 V^A - if A)] 9(-T) = (5.6) 

c 

for all data on an arbitrary compact null C. In particular, we can choose as data 

01s 0*0 = S s (x,x) , (5.7) 

for any nonsingular point x G C. Since x is arbitrary, this implies that either A(x', x) = 0, 
and we get HP directly, or 

V a n a (x) + n a (x)[2V a \n\A(x,x')\- H a (x)] =0 (5.8) 

for all pairs of points with x' E I + (x), x E C nonsingular, and all null vector fields n a 
which are gradients of a function. 
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We will now study the consequences of (5.8), by writing it in a suitable coordinate 
system. In terms of w, (5.8) is of the form 

9 ab V a V b w + X a V a w = , (5.9) 

where X a (x,x') := 2 V a In | A(x, x') \ — H a (x). Choose any point x G C, and a system 
of Riemann normal coordinates based at x , for which g^ u (x) = rj + 0(2), where we 
have introduced the notation 0(n) for a quantity which is of order n in the coordinate 
separation of x from x ; then (5.9) becomes 

rj^d^w + X^d^w = 0(1) . (5.10) 

Let us now consider the function 



w 



= t + yf¥x~ + c^x»x v + 0(3) , (5.11) 



where the c^ v are appropriate coefficients such that g^d^w d v w = 0. Substituting into 
(5.10) we get that, for x l x i ^ 0, 

m — 2 . ri--X' l (x,x')x j 

^= + 2 + X ^ x ) + 3 /— — r = 0{l) . (5.12) 

The first term in the left hand side of (5.12) is 0{— 1), whereas the others are all O(0); 
therefore (5.12) implies m = 2 and 

X\x Q , x') + X x (x , x') hm A = -2 ?TV ' ( 5 - 13 ) 

x^>0 ± \X\ 

where the ± in the limit corresponds to the arbitrariness in the direction of approach of x^ 
to the origin and, with some abuse of notation, we have denoted by x the spatial coordinate 
of the two-dimensional Minkowski spacetime. This means that 

X*(x 0i x') ± X*(x , x') = -2 rT V (5-14) 

must hold simultaneously for both signs, i.e., that X x (x , x') = 0. Since this is true in any 
system of local Minkowskian coordinates, and taking into account the arbitrariness of x , 
we deduce that X a (x, x') = field. 

In conclusion, we have found that CPP implies either HP, or m = 2 and 

H a (x) = 2V a ln\A(x,x')\. (5.15) 

In other words, CPP is always equivalent, for m > 2, to HP; in two spacetime dimensions 
A ^ and HP is always violated, whereas CPP can be satisfied. Therefore, CPP can 
be regarded as a nontrivial generalisation of HP, although it becomes interesting by itself 
only for m = 2. More explicitly, the definitions of CPP and HP differed in that for the 
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former, the support of the data had to be null rather than achronal, and the data had to 
vanish at non-smooth points of C. We now know that this is of no consequence other than 
in some two-dimensional cases, but it is perhaps surprising that such a small restriction 
in the data can have the effect of enlarging, in a nontrivial way, the class of equations to 
which the definition applies. 

We have derived (5.15) above as a necessary condition for the validity of CPP in two 
dimensions, but the following simple argument shows that it is also sufficient. Substituting 
(5.15) into (5.9) the latter reduces to g ah V a V b w = 0, and CPP holds if this equation is 
satisfied by any function w that generates null hyper surfaces. For m = 2 this is true, 
because w is null iff it depends on only one of the coordinates u and v defined by (2.5) 
and (2.6), i.e., iff either w = w^u) or w = w 2 (v). Using the identity (2.2) for / = w, we 
have g ab V a V b w = 0. 

Let us now see what restrictions the condition (5.15) places on the wave equation. 
First of all, since the Riemannian connection is torsion- free, we have Vj a V b ] In | A| =0, and 
(5.15) implies 

V [a H b] =0; (5.16) 

this is an integrability condition that allows us to obtain A by straightforward integration 
of (5.15). In fact, (5.16) implies that H a = 2V a A, for some A; substituting back into 
(5.15), we have 

^■^^W-y <«■»> 

since A(x',x') = 1 by (2.9) and (2.10). Furthermore, it is not difficult to verify, from the 
definition of the Green function, that A(x,x') satisfies the "adjoint wave equation" in x 
(see also (2.8) and (2.11)) 

g ab V a V b A-V a (H a A)+KA = 0. (5.18) 
Substituting (5.15) twice into (5.18), we get 

2 V a H a + g ab H a H b - 4 K = , (5.19) 

hence 

g ab V a V b A + g ab V a AV b A - K = . (5.20) 
Using (5.17), we can replace (5.20) by an equation for A simpler than (5.18), 

g ab V a V b A - KA = , (5.21) 

which can also be directly obtained by using (5.19) into (5.18). Equations (5.16) and (5.20) 
(or (5.19)) are also sufficient conditions for CPP, as one can easily see by considering a new 
field ifj := (j)exp A. From (1.3) and (5.20), we find that ip obeys the equation (7 a6 V a V 6 -i/> = 0, 
which is CPP because m = 2. Since ip and <j> differ only by the nonvanishing factor expA, 
it follows that 4> satisfies CPP as well. 
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In the coordinates (u,v) we have H u = -0" 2 t/, H v = -0" 2 F, K = -fi" 2 W, so 
that (5.16) and (5.19) are equivalent to 

d u U = d v V (5.22) 

and 

d u U + UV-W = 0, (5.23) 
respectively. Similarly, we have U = d v A and V = d u A; (5.20) and (5.21) become 

d 2 uv A + d u Ad v A-W = (5.24) 

and 

8l v A -WA = 0. (5.25) 

These conditions for CPP represent strong constraints on the wave equation; for example, 
when U = V = 0, (5.23) and (5.24) are satisfied only if W = [15]. Notice that when 
(5.16), or equivalently (5.22), is satisfied, the functions a and r defined in (3.6) and (3.7) 
both reduce to A; upon substitution into (3.8)-(3.10), equations (5.22) and (5.23) then 
imply that j = l = 1 and j 1 = l_ x = 0. 

6. Equivalence between CPP and PW in two spacetime dimensions 

In this section we will show that, in two spacetime dimensions, CPP and PW are equivalent 
properties; the proof is very simple if we rely on some results of the previous section. In 
fact, since the two-dimensional wave equation (2.7) satisfies CPP iff <j> can be written as 
4> = exp(— A) ip, with d^ v ip = 0, it follows that (2.7) is CPP iff its general solution has the 
form 

<f>(u, v) = e~ Aiu > v) (r(u) + s(v)) , (6.1) 
where r and s are arbitrary functions. By comparing (6.1) with the = case of (3.3), 

(j>(u,v) = f (u,v) R(u) + g (u,v) S(v) , (6.2) 

we conclude immediately that CPP =^ PW . In order to prove the converse, i.e., that 
PW =^ CPP, substitute the expression (6.2) into (2.7); invoking the arbitrariness of R 
and S, we get: 

; (6.3) 
; (6.4) 

(6.5) 
(6.6) 



dlJo + UdJ + VdJ + Wf = 

d L% + ud u9o + vd v9o + w % = 
dj + Uf = 0; 
d u 9 + Vg = 0. 
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After some manipulations, these equations lead to (5.22) and (5.23), which are sufficient 
conditions for the validity of CPP, as we saw in section 5. We can therefore conclude that, 
in a two-dimensional spacetime, CPP and PW are equivalent. 

An immediate corollary of this result is that not every expression of the type (6.2) is 
the general solution of a two-dimensional wave equation. This is evident by a comparison 
between (6.1) and (6.2), which shows that the amplitudes f and g must be related by 
f (u,v)/g (u,v) = a{u)/[3{v), for some functions a and (3. The same conclusion can be 
reached in a more explicit way by integration of (6.5) and (6.6), to obtain 



with a(u) := f (u, 0) exp A(w, 0) and j3{v) := g (0, v) exp A(0, v). By a generalisation of 
this argument, one could show that not every expression of the type (3.3), with arbitrary 
R and S, is the general solution of a two-dimensional wave equation. 

It is possible to prove that CPP <^ PW without necessarily using the results of the 
previous section. This is obvious for what concerns the implication PW =^> CPP, since 
PW leads us to the same equations (5.22) and (5.23) of section 5, and one can proceed to 
show that they are sufficient conditions for CPP exactly in the same way as we did there. 
The proof that CPP =>- PW is less compact, but we nevertheless present it because of its 
intrinsic interest. 

Consider the null data of compact support given by tp(u) = 5(u — -u ), with u > 0, 
and ip(v) = 0; then CPP implies that <J)(u,v) must be concentrated on the line u = u . 
This means that, for any fixed v, (f)(u, v) must be a distribution in u that, applied to a 
test function F(u, v), produces a number depending only on the behaviour of F at the 
point (u ,v), i.e., on fW(u ,v), with i > 0. Since distributions are linear functionals, the 
combination of the various F^>(u , v) must be linear, so that we can write 



for some functions f i . The sum on the right hand side of (6.9) actually consists only of 
a finite number of terms; this follows from a rigorous result of distribution theory [23], 
but can be heuristically justified as follows. Since the action of as a functional must be 
defined on any test function F, the value of the series 



must be finite for each F. Suppose that the sum in (6.9) and (6.10) is infinite, that is, that 
there exist infinitely many f i ^ 0. Then, taking F such that F^ l >(u , v) = /^Wq, i>) -1 when 




(6.7) 



(6.8) 




(6.9) 




(6.10) 
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f i 7^ 0, and arbitrary otherwise, the expression (6.10) becomes infinite, which contradicts 
the hypothesis of regularity of <fi. Hence, we must have 



(f)(u,v) = J2fi( u O' v ) s(l) ( u - u o) > 



(6.11) 



i=0 



with N u finite. Similarly, if <p(u) = and ip(v) = S(v — v ), we have 



N v 



(j)(u,v) = ^^(m,v )5 w (v -v ) 



(6.12) 



i=0 



for some finite N v and some functions g i . 

From (6.11) and (6.12) it already follows that (2.7) is PW^, with N = max{iV u , N v }, 
as can be seen just by using linearity and representing arbitrary data <p and ip as super- 
positions of <5- like data, 



(p(u) = J du <p(u ) 8(u - u Q ) , 
^( v ) = J dv ip(v )S{v -v ) , 



(6.13) 
(6.14) 



but without making any other use of the differential equation for <p. We shall now show, 
however, that (2.7) allows us to set N = 0, i.e., to conclude that CPP =>- PW . For this 
purpose, let us first consider again data <p(u) = S(u — u ), ip(v) = 0. Substituting the 
solution (6.11) into (2.7) we get 

[d v f Nu (u ,v) + U(u,v)f Nu (u ,v)] 8^ +1 \u-u ) + 

N u 

+ [V(u, v) dj> , v) + W{u, w )/ 4 (« , v) + 

i=l 

+ d v f i _ 1 (u ,v) + U(u,v)f i _ 1 (u ,v) 5 {l) (u-u ) + 

V(u , v) dj (u , v) + W(u Q , v) f Q {u , v) \ 8{u -u ) = 0. 



+ 



(6.15) 



Smoothing (6.15) with a test function, and taking into account the arbitrariness of the 
latter, we obtain the following set of equations: 



d vfNAu ,v) + U(u ,v)f N (u ,v) = ; 



(6.16) 



i=k 
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+ di- k U(u ,v)f i _ 1 (u ,v) +(-1) 



N u + 1 / iv ii 1" \ f) N u -k+l 




U(u ,v) f Nu (u ,v) = ; 



N u > k > 1 ; (6.17) 




(6.18) 



From (6.16) and the fact that f N (u , 0) = 0, which follows from using our data in (6.11), 
we obtain that f N (u , v) = for all v. Using this in (6.17) and repeating the procedure, 
we obtain f^u^.v) = 0, Mi > 1, so that only f (u ,v) can be nonzero. An analogous 
argument, using data <p(u) = 0, ip(v) = S(v — v ), shows that the only nonvanishing 
coefficient in (6.12) can be g (u, v ). Therefore, the general solution, obtained from (6.11) 
and (6.12) using arbitrary initial data as in (6.13) and (6.14), is 



from which the validity of PW follows. Furthermore, it is not difficult to see that (6.17) 
for k = 1 reduces to (6.5) in u = u , and that substituting this relation into (6.18) we 
obtain (5.23). By a completely symmetric procedure we can recover (6.6). 

7. Equivalence between PWjy and KN^ 

As we saw in section 3, all KN^ wave equations are PW^; it has been conjectured that 
the converse is also true, namely that all ~PW N equations can in fact be solved exactly 
by the Kundt-Newman method in steps [20,24]. To show that this is indeed the case, 
suppose we have a 2-dimensional wave equation whose general solution is the progressing 
wave (3.3). By suitable transformations, we can always write the wave equation in its th 
v- or -u-normal form, (3.4) or (3.5); let us concentrate on the ^-normal form. 

Applying the differential operator in (3.4) to the progressing wave (3.3), and imposing 
that the coefficients of derivatives of the arbitrary functions R(u) and S(v) of different 
orders vanish separately in the resulting expression, we have: 



(f)(u, v) = f (u, v) <p(u) + g (u, v) %j)(v) , 



(6.19) 



d v (j d u f )-j 1 f = , 

d vUo d u fi) ~ 3i fi + d v Uo fi-i) = . 

d v(jof N ) = . 



1 < z < N , 



(7.1) 
(7.2) 
(7.3) 
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and 



d vtio d u9o) ~ 3i 9 = • 

d vUo d u 9i) - k 9t + Jo d u9i-i = , 1 < i < N , 
d u 9 N = , 



(7.4) 
(7.5) 
(7.6) 



Let us consider the group (7.4)-(7.6). Equations (7.4) and (7.5) can be rewritten more 
conveniently as 



Ji 



d v [-du9o) +d v ln\ 3l \ J -°d u g -g = 



Ji 



(7.7) 



and 



d v [ f du9i J + d v In \ h | f d u9i - 9i + f d u9l _ x = , 1 < i < N , (7.1 



respectively. By repeated differentiation of (7.8) with respect to u, one finds that, for any 
h > 1 and k = 0, .. ., N - 1, 

d v V h u g N _ k + d v In \j h \ V h u g N _ k - V^g^, + V h u g N _ k _ x = , (7.9) 
where the differential operator X>£, containing k derivatives d u , is defined as 



Jfc \Jk-l V ^2 VJl 



(7.10) 



for any sufficiently differentiable function F. We shall prove (7.9) by induction over h. 
Assuming that it holds for some h, and differentiating it with respect to u, we have 

dl v V h u g N _ k + d v In \j h \ d u V h u g N _ k - J -f± V h u g N _ k + d u V h u g N _ k _ 1 = , (7.11) 

Jh 

where we have used (3.11). Rewriting the first term as 



i^-^v^;.</.y a ) 

Jh Jh+1 



Jh+1 



Jh 



d u V h u g N _ k + ^ d v V h u +1 g N _ k , (7.12) 
Jh 



we see that (7.11) reduces to the form taken by (7.9) when h — > h + 1. Moreover, for h = 1 
(7.9) becomes 

d v d u9N _ k ) + d v In \ h \ J f d u g N _ k - g N _ k + 3 f d u g N _ k _ 1 = , (7.13) 
\Ji / Ji Ji 

which is just (7.8) for i = N — k, hence certainly true. This completes the proof of (7.9). 

19 



In the proof that the substitution sequence terminates, we need (7.9) only to derive 
the fact that, for < k < N, 

g N = v k u g N - k . (7.14) 

The proof of (7.14) is also by induction, this time over k. Assuming that the equation 
holds for some k < N — 1, and writing (7.9) for the particular case h = k + 1, we have 

d v [ir- 9 u 9n) + d v In |j fc+ il -I*- d u g N - g N + V k + 1 g N _ k _ 1 = . (7.15) 
\Jk+i J Jk+i 

Now, (7.6) allows us to obtain 

g N =V k u + 1 g N _ k _ 1 , (7.16) 

i.e., (7.14) for k — > k + 1. Since for k = (7.14) is trivially true (and even for k = 1 it can 
be easily obtained by substituting (7.6) into (7.8) for i = N), its validity for < k < N is 
established. In particular, we shall be interested in the case k = N, for which 

9n = ^u9o- (7-17) 

The equation analogous to (7.9) with k = N is found by repeated differentiation of 
(7.7) with respect to u; we have, for i > 1, 

d v V\ g + d v In VI g - V^ 1 g = , (7.18) 

whose proof is perfectly analogous to that of (7.9). For i = N, (7.18) becomes, using 
(7.17), 

d v g N + d v In \j N \ g N - V^g, = . (7.19) 

The final step in the proof consists in taking a further derivative with respect to u of 
(7.19), which gives, by (3.11), (7.6), and (7.17), 

^^ = 0, (7.20) 

Jn 

i.e., j N+1 = 0. We have thus shown that the substitution sequence of a PW^ wave equation 
is upper terminating in iV steps. To prove that it is also lower terminating in iV steps, we 
could manipulate (7.1)-(7.3) to show that j_ N _ x = oo, i.e., l_ N _ x = by (3.17). However, 
it is much easier to notice that, starting with the u- normal form of the wave equation, we 
can repeat the proof above in a completely symmetric way to get directly l_ N _ x = 0. 
Hence, PW W =^ KN N and, since we know already that the converse is also true, we can 
conclude that PW^ and KN^ are equivalent properties. 
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8. Conclusions and open questions 

The results presented in this paper clarify the relationships between several properties 
of wave equations related to the absence of tails in their solutions, some of which were 
obvious, while some others were never explicitly analysed in the literature, at least to 
our knowledge. We have seen that the Huygens principle (HP) and the tail-free property 
(TF) are satisfied by the same equations, that the characteristic propagation property 
(CPP) is more general only in that it is satisfied, in addition, by special two-dimensional 
equations, and that the progressing wave propagation property (PW N ) and the solvability 
by the Kundt-Newman method (KNjy) for two-dimensional wave equations are equivalent. 
The two latter properties were not defined in more than two dimensions, and they are 
not given in geometrical terms, despite the motivation we gave for introducing PW N in 
section 3. However, we also showed that in two dimensions PW is equivalent to CPP, 
and acquires thus a geometrical meaning. It would be therefore interesting to investigate 
possible geometrical aspects of PW W equations for higher iV, and whether the definition of 
PW^ can be meaningfully extended to higher dimensions; we will comment here on these 
issues. 

The result of section 6, that in 1 + 1 dimensions CPP is equivalent to PW , has the 
obvious corollary that no ~PW N equation with N > can satisfy CPP. This consequence 
is at first surprising, because from the general form (3.3) of a progressing wave, one is 
tempted to think that, in two dimensions at least, all PW W wave equations satisfy the 
CPP. For, by choosing R(u) and S(v) to be of compact support, their derivatives will 
also be of compact support, and all of will be made of pieces sandwiched between null 
coordinate lines, so it will have no tails. The problem with this argument is that, although 
solutions with no tails are indeed obtained when R and S are of compact support, the 
latter functions are not themselves the null data which are involved in the definition of 
CPP; rather, the data are 

N N 

i=0 i=0 
N N 

i=0 i=0 

and the R(u) and S(v) that correspond to generic <p and ip of compact support, will in 
general involve integrals of <p and ip if N > 0, and will thus be of non-compact support. 

We can clarify this point by considering a typical example of PW^ equations, ob- 
tained by separating the angular variables in the massless Klein-Gordon equation □ $ = 
in 4-dimensional Minkowski spacetime. Expanding $ in spherical harmonics as $(x) = 
Yliim Xim(t> r ) Yim{Q-i </?), and rescaling the radial components to 4> lm := r Xi m -> we have that 
the 4> lm satisfy the /-dependent equation 

<t>lm = , (8-3) 



d? - dl 



+ 



1(1 + 1) 



Im 



8 2 , 1(1 + 1) 

uv (v- u) 2 
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where the null coordinates u and v are defined as in (2.5) and (2.6), but now with x replaced 
by r; the two-dimensional equation (8.3) is PW ; and can be solved exactly [17,20]. Since 
(5.23) is a necessary condition for CPP, and in the case of (8.3) we have U = V = 0, it 
follows that CPP holds only if W = 0; this is clearly false for I > 1. Let us see explicitly 
the failure of CPP in the simple case 1 = 1. Equation (8.3) becomes 



d 2 +— 

uv (v- u) 2 



4>(u, v) = 



(8.4) 



which is PW X ; its general solution is 

2R{u) 



(f>(u,v) 



v — u 



+ R'(u)-^- + S'(v) 



v — u 



(8.5) 



We now look for the particular solution generated by null data with support at one point, 
e.g., 



ip{u) = 8{u - u Q ) , 
ip(v) = . 



(8.6) 
(8.7) 



These data correspond to R(u) and S(v) satisfying the differential equations 



R'(u) - - R(u) + - S(0) + S'(0) = S(u - u ) 

%Jb Xii 



(8.8) 



and 



S'(v) - - S(v) + - R(0) + R'(0) = , 

v v 



whose general solutions are, respectively, 

,2 



U 



R(u) = —6(u- u ) + au 2 + S'(0) u + S(0) 



u 







(8.9) 



(8.10) 



and 



S(v) = b v 2 + R'(0)v + R(0) 



(8.11) 



with a and b arbitrary constants. Substituting (8.10) and (8.11) into (8.5), and using the 
relations R(0) = S(0) and R'(0) = S'(0) that follow, e.g., from (8.11), we have 



, / \ c/ , 2 uv 

<p(u, v) = d(u — u )-\ — 2 



Uq v — u 



6{u - u ) + c] , 



(8.12) 



where c := (a — b) Uq can be prescribed arbitrarily. It is clear from this expression that 
the only line v = const on which is of compact support is v = 0; the choice c = gives 



22 



a retarded wave, c = —1 an advanced one, and any other value gives a solution which is 
non-zero almost everywhere on Ail 

Turning now to the question of whether it is meaningful to define the class of wave 
equations whose general solution is a finite sum of progressing waves, in dimension m > 2, 
the simplest candidate for such an equation would again be the example □<& = in 4- 
dimensional Minkowski spacetime, which satisfies HP. Using the expansion in spherical 
harmonics given above, and (8.3), whose general solution is of the type (3.3) with N = I 
and all the functions carrying labels I and m, we can write 

+ 00 11/ \ 

= E E E ( r filmic r) Y lm (0, <p) R®(t -r) + rg ilm (t, r) Y lm (e, <p) S®(t + r) , 

1=0 i=0 m=-l ^ ' 

(8.13) 

which contains derivatives up to arbitrarily high order of an infinite number of free functions 
R lm and S lm , and is therefore not of the form (3.3). This shows that, at least in four 
dimensions, the definition of PW^ as it stands is empty for the case of spherical wavefronts, 
and makes it plausible that it is ill-posed in general. 

It is still useful to define a "reduced" PW^ property for an m-dimensional equation, 
in the sense that all two-dimensional equations, obtained by separating out appropriately 
chosen coordinates, may have progressing wave general solutions, as in the case of (8.3). A 
more general example, in four dimensions, would be obtained with equations of the type 

f{y, z) K t x + g(t, x) H y ^ <f>(t, x,y,z) = 0, (8.14) 

where / and g are nonvanishing functions and K t , H are suitable differential op- 
erators acting on the variables used as subscripts. With the separation of variables 
(f)(t,x,y,z) = tp(t,x) cp(y, z), we see that (8.14) satisfies the reduced progressing wave 
propagation property if, for each value of a admitted by the equation 



H y z <p(y, z) + a f(y, z) <p(y, z) = , (8.15) 

the reduced wave equation 

K t>x if>(t, x)-a g(t, x) x) = (8.16) 



is PWjy for some N(a). The high degree of dependence of the reduction procedure — hence 
of the properties of the reduced equation — from specific, noncovariant structure is evident. 
It is thus very difficult to make general statements about this property and its possible 
relationships with other ones we have discussed in this paper, independently of the explicit 
form of the wave equation and the choice of coordinates singled out for separation (i.e., of 
wave front). 

It would be interesting to know whether a meaningful, truly m-dimensional generali- 
sation of the PW^ property can be given, for another reason as well. The expressions for 
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the advanced or retarded Green function of = in m-dimensional Minkowski space, 
which for even m > 4 are of the form [5,7] 

N 

G m (t',x';t,x) = 5> mi (|x' - x|) a«(f - * ± | X ' - x|) , (8.17) 

with N = (m — 4)/2, remind one of definition (3.2), and suggest a possible connection 
between HP and progressing waves, this time of order iV > 0, in dimension m > 2. So far 
we have not been able to find such a relationship, however. For now, ~PW N equations are 
left with no simple interpretation in terms of tails, and no true generalisation to higher 
dimensions, but their usefulness derives from the fact that they are exactly solvable by the 
method of Kundt and Newman [20] , as we saw in section 7, and from their relationship to 
To da lattices [25]. 

Taking into account all of these remarks, there is little doubt that, as far as studies 
of tails are concerned, the definition that should be preferred is, for its generality, the one 
based on HP or, equivalently, on TF. With a notion of wave tails that is now completely 
clear and unambiguous, one can conduct further research in several directions. Here is a 
list of some problems that one may address: 

(i) Which are the conditions on g ab , H a , and K that make (1.3) tail-free? Because of the 
geometrical interpretation of g ab as a metric in A4, it is natural to restate this question 
by asking in which spacetimes (i.e., for which class of g ab ) is the wave equation tail- 
free. This is sometimes called the Hadamard problem, and has received some attention 
over the years [6,26,27]. 

(ii) What happens for vector and tensor fields? These cases include the propagation 
of electromagnetic and gravitational waves [15,26,28] and are more realistic — though 
somewhat more complicated — than the scalar field one. 

(iii) Which are the physical consequences and effects of wave tails? The situation that 
has been investigated most in detail is that of radiation from compact objects [2,3,8]. 
The corresponding problem in a cosmo logical background [11] has apparently been ne- 
glected, probably because of the belief that any observable effect should be extremely 
small. However, in this case curvature never drops off and scattering continues forever, 
so it might convey a relevant amount of radiation in the interior of the light cone [9] . 
This brings immediately out a new problem, namely how much radiation goes into 
the tail, which leads us to the next question. 

(iv) How to quantify tails? A very natural and intuitive way would be to calculate a reflec- 
tion coefficient; this is possible when backscattering is localised and there are regions 
in which the field is free (even asymptotically, or after a suitable coordinates trans- 
formation has been performed [8]). In a cosmological context, however, the property 
that makes the phenomenon potentially interesting — i.e., the fact that backscattering 
takes place always and everywhere — at the same time prevents us from defining purely 
ingoing and outgoing solutions of the wave equation, and hence from computing re- 
flection and transmission coefficients [9]. It is necessary to find alternative ways to 
quantify tails, perhaps based on the ratio between their energy content and the total 
energy of the field. 



24 



Acknowledgements 

This work was supported by the Commission of the European Communities under the 
DG-XII contract no. CIl*-0540-M(TT) and the DG-III contract no. ECRU002, and by 
the Instituts Internationaux de Physique et de Chimie Solvay. 

References 

[1] Huygens C 1690 Traite de la lumiere (Leiden: Van der Aa) 

[2] Blanchet L and Damour T 1992 Hereditary effects in gravitational radiation Phys. Rev. D 46 4304-19; 
Blanchet L and Schafer G 1993 Gravitational wave tails and binary star systems Class. Quantum Grav. 
10 2699-721 

[3] Wiseman A G 1993 Coalescing binary systems of compact objects to (post) 5 / 2 -Newtonian order. IV. 

The gravitational wave tail Phys. Rev. D 48 4757-70 
[4] Hadamard J 1952 Lectures on Cauchy's Problem in Linear Partial Differential Equations (New York: 

Dover) 

[5] Courant R and Hilbert D 1962 Methods of Mathematical Physics (New York: Wiley) vol 2 
[6] Friedlander F G 1947 Simple progressive solutions of the wave equation Proc. Cambridge Philos. Soc. 
43 360-73 

[7] Soodak H and Tiersten M S 1993 Wakes and waves in N dimensions Am. J. Phys. 61 395-401 
[8] Price R H 1972 Nonspherical perturbations of relativistic gravitational collapse. I. Scalar and gravita- 
tional perturbations Phys. Rev. D 5 2419-38; 

Misner C W, Thorne K S and Wheeler J A 1973 Gravitation (San Francisco: Freeman) 
[9] Faraoni V and Sonego S 1992 On the tail problem in cosmology Phys. Lett. 170A 413—20 
[10] Visser M 1993 Acoustic propagation in fluids: an unexpected example of Lorentzian geometry (preprint 
gr-qc/9311028| ) 

[11] Ellis G F R and Sciama D W 1972 Global and non-global problems in cosmology, General Relativity, 
Papers in Honour of J. L. Synge ed L O'Raifeartaigh (Oxford: Clarendon) 35-59 

[12] DeWitt B S and Brehme R W 1960 Radiation damping in a gravitational field Ann. Phys. (N.Y.) 9 
220-59 

[13] Fox R et al 1969 Do faster-than-light group velocities imply violation of causality? Nature 223 597; 

1970 Faster-than-light group velocities and causality violation Proc. R. Soc. A 316, 515-24; 

Bers A et al 1971 The impossibility of free tachyons Relativity and Gravitation C G Kuper and A 

Peres (New York: Gordon and Breach) 
[14] Wald R M 1984 General Relativity (Chicago: University of Chicago Press) 

[15] Sonego S and Faraoni V 1992 Huygens' principle and characteristic propagation property for waves in 

curved space-times J. Math. Phys. 33 625-32 
[16] Friedlander F G 1975 The Wave Equation on a Curved Space-time (Cambridge: Cambridge University 

Press) 

[17] Couch W E and Torrence R J 1986 A class of wave equations with progressive wave solutions of finite 

order Phys. Lett. 117A 270-4 
[18] Ward R S 1987 Progressing waves in flat spacetime and in plane-wave solutions Class. Quantum Grav. 

4 775-8 



25 



[19] Torrence R J and Couch W E 1988 Progressing waves on spherical spacetimes Gen. Rel. Grav. 20 
343-58 

[20] Kundt W and Newman E T 1968 Hyperbolic differential equations in two dimensions J. Math. Phys. 
9 2193-210 

[21] Torrence R J and Couch W E 1985 Transparency of de Sitter and anti-de Sitter spacetimes to multipole 

fields Class. Quantum Grav. 2 545-53 
[22] Bombelli L, Couch W E and Torrence R J 1991 Wake-free waves in one and three dimensions J. Math. 

Phys. 32 106-8 

[23] Richtmyer R D 1981 Principles of Advanced Mathematical Physics (New York: Springer- Verlag) vol 
I, P 51 

[24] Torrence R J 1990 Self-adjoint acoustic equations with progressing wave solutions J. Phys. A: Math. 
Gen. 23 4107-15 

[25] Torrence R J 1987 Linear wave equations as motions on a Toda lattice J. Phys. A: Math. Gen. 20 
91-102; 

Bombelli L, Couch W E and Torrence R J 1992 Solvable systems of wave equations and non-Abelian 
Toda lattices J. Phys. A: Math. Gen. 25 1309-27 
[26] Kiinzle H P 1968 Maxwell fields satisfying Huygens's principle Proc. Camb. Phil. Soc. 64 779-85; 

Carminati J and McLenaghan R G 1986 An explicit determination of the Petrov type N space-times 
on which the conformally invariant scalar wave equation satisfies Huygens' principle Ann. Inst. Henri 
Poincare, Phys. Theor. 44 115-53 

1987 An explicit determination of the spacetimes on which the conformally invariant scalar wave 

equation satisfies Huygens' principle. Part II: Petrov type D space-times Ann. Inst. Henri Poincare, 
Phys. Theor. 47 337-54 

1988 An explicit determination of the spacetimes on which the conformally invariant scalar wave 

equation satisfies Huygens' principle. Part III: Petrov type III space-times Ann. Inst. Henri Poincare, 
Phys. Theor. 48 77-96 

[27] McLenaghan R G 1969 An explicit determination of the empty space-times on which the wave equation 
satisfies Huygens' principle Proc. Camb. Phil. Soc. 65 139-55 

[28] Noonan T W 1989a Huygens's principle for the electromagnetic vector potential in Riemannian space- 
times Ap. J. 341 786-95 

1989b Huygens's principle for the wave equation for second-rank tensor fields Ap. J. 343 849-52; 

Wiinsch V 1990 Cauchy's problem and Huygens' principle for the linearized Einstein field equations 
Gen. Rel. Grav. 22 843-62; 

Caldwell R R 1993 Green's functions for gravitational waves in FRW spacetimes Phys. Rev. D 48 
4688-92 



26 



